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We derive expressions for the gluon production cross sections in the single diffractive 
C^l ' proton-nucleus scattering and DIS processes in a quasi-classical approximation. The result- 

ing cross sections include the effects of all multiple rescatterings in the classical background 
>D . field of the target proton or nucleus, which remains intact after the scattering. We also 

write down an expression for the inclusive gluon production cross section in DIS in the 
quasi-classical approximation. 
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In the recent years it has become clear that partonic saturation may play an important and possibly vital 
role in high energy hadronic and nuclear collisions. Partonic saturation, associated with slowing down of the 
growth of structure functions with energy, is characterized by high density of quarks and gluons and strong 
gluonic fields in the hadronic and nuclear wave functions [1-4]. The transition to the saturation region is 
described by the saturation scale Q 2 s {s), which is an increasing function of the center of mass energy s and 
for high enough s could lie in the perturbative QCD region [2,5,6]. It has been suggested that most of the 
quarks and gluons in the small- a; hadronic wave function have the transverse momentum of the order of Q s 
[7], making perturbative calculation possible for most observables in high energy hadronic collisions. 

Calculations in the saturation region consist of two stages: classical and quantum. As it was conjectured 
in [2] the strong gluonic fields in the small- a; tail of a hadronic or nuclear wave function could be very well 
approximated by the classical solution of the Yang-Mills equations of motion. Even though this conjecture 
was originally formulated with the valence quarks as sources of color charge generating the gluon field [2] 
it could be generalized to include all the higher momentum components of the nuclear wave function in 
the source current [8]. The small- a; gluons have a large coherence length of the order of Xjlm^x allowing 
them, for small enough x, to coherently interact with all the nucleons in the nucleus [9]. The classical 
gluonic field of a nucleus, known as the non-Abelian Weizsacker- Williams field, was shown to include all the 
multiple rescatterings of the gluons with the sources of color charge [3,4]. In the valence quark formulation 
of the problem for the nucleus where the multiple rescatterings are mediated by two gluon exchanges this 
parametrically corresponds to resumming all powers of the parameter a 2 A 1 ' 3 , with A the atomic number 
of the nucleus [3] . In terms of dynamical variables resummation of this parameter is equivalent to summing 
up powers of Q 2 s /k\, where k± is the transverse momentum of the gluon. The physical effect of multiple 
rescatterings is to push the gluon's transverse momentum from some almost non-perturbative scale up to 
k±_ ~ Q s thus suppressing the infrared region and the phenomena associated with it [10]. 

The second stage of the saturation calculations corresponds to inclusion of quantum corrections. Usually 
the quantum corrections are taken in the leading logarithmic approximation and bring in powers of a s In 1/x, 
though in principle subleading logarithmic contributions should be considered, since while being paramet- 
rically smaller they could still be of numerical importance [11]. In the traditional language the leading 



logarithmic corrections correspond to resummation of multiple BFKL pomeron exchanges [12,1]. They can 
also be represented as a series of emissions of the non-Abelian Wcizsacker- Williams field, where after each 
interaction the produced field is incorporated in the hard source for production of softer field at the next 
step [8]. 

Several observables have been calculated using the saturation approach. The one which probably received 
most attention is the total cross section of deep inelastic scattering (DIS), since it is related to quark and 
gluon distribution functions. In the "classical" approximation it has been first calculated for QCD in [13], 
yielding a Glauber-type formula. Then an extensive effort was put into inclusion of the quantum corrections 
in the total DIS cross section and the gluon distribution function resulting in the nonlinear equation of [5,6]. 
The equation could be derived in high energy effective theory of [6] or using the dipolc model of [14]. Later 
the nonlinear equation has been reproduced and solved by approximate and numerical methods in [15]. 

Another important observable which could be calculated in the saturation framework is the inclusive single 
gluon production cross section in DIS, proton-nucleus (pA) and nucleus- nucleus (AA) collisions. The gluon 
production cross section could be compared with the data on minijet and pion production at mid-rapidity. 
In the first stage of the saturation calculation the problem of finding the inclusive single gluon cross section 
is equivalent to finding the classical field produced by two colliding nuclei or hadrons which was found at the 
lowest order in [16]. The problem of inclusive gluon production in pA including all orders of rescatterings in 
the nucleus, i.e., resumming all powers of a 2 s A 1 ' 3 , has been solved in [17] and the result was later reproduced 
in [18,19]. At the same time the question of inclusive gluon production in DIS has not been studied yet. The 
calculation there is a little different from pA and will be discussed later in this paper. In AA collisions the 
problem is more complicated than the single nucleus case of DIS since now one has to include rescatterings 
of the produced gluon in the fields of nucleons in both nuclei. The gluon production problem in AA is of 
utmost importance for initial conditions for quark-gluon plasma (QGP) formation in heavy ion collisions. An 
important recent progress has been made in [20-23] on understanding the distribution of produced gluons in 
AA and in [24] on exploring the possible thermalization of gluons in the saturation model. At the moment 
nobody found a way to include the quantum corrections to the classical inclusive gluon production cross 
section neither in DIS nor in pA or AA and this important question still remains open. 

The issue of multiple rescatterings in a medium produced in the heavy ion collisions is related to the 
problem of energy loss of produced particles [25]. Recently there has been developed a vigorous activity 
in the field of energy loss [26-28] due to sensitivity of this observable to the possible QGP formation in 
heavy ion collisions. Multiple rescatterings in the cold nuclear medium [26] are of similar nature to multiple 
rescatterings of the saturation calculations [17] and the interplay of two notions could be very important for 
our understanding of nuclear collisions. 

Exclusive observables, such as diffractive cross section and particle production are also very important 
for our understanding of high energy scattering. The total diffractive cross section and the corresponding 
structure function have been found for DIS in the quasi-classical approximation in [29-31]. A one-gluon 
correction to the total diffractive cross section has been found in [31]. The second stage of the saturation 
calculations has already been completed for diffraction: in [32] an equation has been written resumming all 
leading logarithmic corrections (multiple pomeron exchanges) to the total diffractive cross section in DIS. 
Exclusive vector meson production including all the multiple rescatterings at the quasi-classical level in DIS 
has been studied in [33,31]. 

Nevertheless the question of calculating the diffractive gluon production cross section in the quasi-classical 
approximation has not been addressed yet in the literature. The problem is very important for the proton- 
nucleus collision experiments at RHIC and for DIS at HERA and we are going to study it in this paper. We 
begin in Sect. IIA by considering the case of a quarkonium scattering on a target nucleus. We are interested 
in the gluon production process where the nucleus remains intact with no restrictions imposed on possible 
final states of the qq pair. The resulting gluon production cross section is given by Eqs. (11) and (17) and 
includes the effects of all multiple rescatterings in the background fields of the nucleons in the nucleus in 
the quasi-classical approximation (equivalent to resummation of powers of a 2 A 1 ' 3 ). We then continue in 
Sect. IIB by applying the developed techniques to the case of proton-nucleus collisions, where the proton 
is approximated as a color single state of three valence quarks. The calculation for the pA case is a little 
more involved than the calculation for quarkonium-nucleus scattering. The answer in a more compact form 
is given by Eqs. (29), (33), (34), (35). We argue that the presence of intrinsic sea quarks and gluons in 
the projectile's wave function may tremendously complicate the calculation of diffractive production cross 
sections, while posing no threat to inclusive cross sections. In Sect. Ill we generalize the result of Sect. IIA 



to the case of deep inelastic scattering. We discuss the differences between inclusive gluon production in DIS 
and pA, which are rooted in the model of the proton one uses. We conclude the paper by writing down an 
expression for inclusive gluon production in DIS given by Eq. (43). 

II. HADRON-NUCLEUS COLLISIONS 

In this section we will first derive an expression for diffractive gluon production cross section in the 
quarkonium-nucleus scattering. We will continue by generalizing the result to the case of diffractive proton- 
nucleus scattering. 

A. Quarkonium— Nucleus Collisions 

Let us start by considering a scattering of a qq (quarkonium) state on a target nucleus. We want to calculate 
diffractive gluon production cross section. That is we are interested in the hadron-nucleus scattering processes 
where a soft (small- x) gluon is produced in the central rapidity region while the nucleus remains intact in the 
color neutral state. The diffractive gluon production cross section could be written as a convolution of the 
onium's wave function squared with the diffractive gluon production cross section for the qq pair scattering 
on a nucleus 

-if^T = J ^M^-ir_,«)\^ y{L) , (1) 

where r is the transverse separation of the quarks in the onium and a is the light cone momentum fraction 
carried by the quark. 

quarkonium 





nucleus 



B 

FIG. 1. Diagrams contributing to the diffractive gluon production cross section in the quarkonium-nucleus scat- 
tering. 

We are going to perform our calculations of da^f , , / d 2 k dy in the time-ordered light cone perturbation 
theory [34] working in the light cone gauge of the projectile hadron A+ = [17]. Similarly to [17,19,21,30,31] 
we distinguish the cases when the gluon is present in the quarkonium's wave function before the collision 
and when it is emitted after the collision. The diagrams contributing to the diffractive gluon production are 
shown in Fig. 1. We assume that the quark and the antiquark in the quarkonium are in the color singlet 
state, that is there is no soft non-perturbative partons which may carry some of the color in the quarkonium 
wave function making the qq pair non-singlct in color. The graph in Fig. 1A represents the case when the 
qq state develops a gluon fluctuation long before hitting the target. Subsequently quark, antiquark and the 
gluon rescatter in the target nucleus leaving it intact since we are interested in diffractive processes [29,30]. 
The exchanged gluons coming from the nucleons in the nucleus can connect to either quark, antiquark or 
gluon lines and we have to sum over all of these possibilities. This is represented in Fig. 1 by not connecting 
gluon lines specifically to any of the projectile parton lines. The diagram in Fig. IB corresponds to the case 
when only qq pair rescatters in the nucleus and the gluon is emitted after the scattering. Again both the 



quark and the antiquark lines interact with the nucleus which is represented by the gluon lines not connecting 
to any specific quark line. Finally in the diagrams of Fig. 1 the gluon is shown to be emitted off one of the 
quark lines only. To obtain the answer one of course has to sum over all possible gluon emissions off the 
quark and antiquark lines in the amplitude and in the complex conjugate amplitude. 

To calculate the diagrams in Fig. 1 one needs to know forward "propagators" of the qq and qqG states 
through the nucleus. The propagator of the quark- antiquark state through the nucleus is well known and is 
given by [13,17,21,30,31] 

^(-^{x-yfQ^ (2 ) 

where x and y are the transverse coordinates of the quark and the antiquark correspondingly and we have 
used a shorthand notation [26,17,21] 



(x-y_) 2 Q 2 s = (3L-V) 2 *-^0 * pxG(x,l/\x-yf), (3) 

with p = yl/[(4/3)7r_R 3 ] the density of the atomic number A and b the impact parameter of the qq pair. 
Here the nucleus is assumed to be a sphere of radius R, but the result could be easily generalized to other 
geometries. The gluon distribution function in Eq. (3) at the two gluon level is given by [26] 

xG(x,l/x 2 ) = ^£ m -L (4) 
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with p, some infrared cutoff. In principle the saturation scale Q s has to be found from the implicit equation 
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QsW = n 2 ~1 pxG(x,Q s (b)). (5) 

Since the gluon distribution of Eq. (4) is a slowly (logarithmically) varying function of transverse distance 
in certain problems the Q s dependence on the right hand side of Eq. (5) could be neglected turning Eq. (5) 
from equation into an equality. 

The "propagator" of the qqG state through the nucleus is a little harder to calculate. To do this let us 
consider the interaction between the qqG state and the first nucleon in the nucleus that exchanges gluons 
with it. The diagrams wc need to sum are depicted in Fig. 2. To write down the contributions of these 
diagrams let us first introduce the gluon-nucleon interaction potential as a Fourier transform into coordinate 
space of the normalized gluon-nucleon scattering amplitude — jSi [26,17] 

With the help of Eq. (6) one can easily see that the contributions of diagrams in Fig. 2 are 

A = ~Wj® (7a) 

B = -§- c V(0) (7b) 

d = §-vU-v) (7d) 
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(71) 



where we take a trace in the color space of the nucleons and average over the colors of the nuclcon, which 
yields a factor of 1/2N C . 
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FIG. 2. Diagrams included in a diffractive scattering of the qqG state on a single nucleon in the target nucleus. 



Summing up the diagrams of Fig. 2 we obtain 

A + B + ... + F = -\[V{Q)~V(x-z)] - ±[V(0)-V(y-z)] + ^[V(0)-V(x-y)]. (8) 

The calculation which led to Eq. (8) never assumed summation over the colors of the quark-antiquark pair 
and the gluon. Eqs. (7) explicitly show us that diffractive scattering on a nucleon does not change the color 
structure of the incoming qqG state. This allows us to easily see that the full answer for the "propagator" 
of the qqG state through the nucleus could be obtained b y just ex ponentiating the interaction with a single 
nucleon. Multiplying the right hand side of Eq. (8) by 2\/R 2 — b 2 pa to take into account the density of the 
nucleons in the nucleus [17] yields 

■ + F) = -\{*- z?Ql-\{iL- z?Ql + ^{x-y) 2 Ql (9) 



-P(x,y,z) = 2^R 2 -b 2 pa(A + B 



Exponentiating Eq. (9) we obtain the propagator of the qqG state through the nucleus 



exp(-P(x,j/,z)) = exp(--(x-z) 2 Q 2 s - - 



te-tfoi + Mpte-tdoi,- 



(10) 



which agrees with the result derived by Kopeliovich et al in [18] and by Kovner and Wiedemann in [31]. 

Now it is straightforward to write down the cross section for diffractive gluon production. The cross 
section is given by the square of the amplitude depicted in Fig. 1 in momentum space. Fixing the overall 
normalization yields 



d ° q d!ff 

d 2 kdy it 2 (2tt) 2 



as ° F ' ' d 2 bd 2 Zl d 2 z 2 e-^-^ ' " ' 
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where Zi and z 2 are the transverse coordinates of the gluon in the amplitude and in the complex conjugate 
amplitude correspondingly. Eq. (11) gives us the cross section of diffractive gluon production in quarkonium- 
nucleus collisions. In general the integrals over z\ and zi in Eq. (11) are very hard to do analytically and they 
should probably be done numerically. However, for not very large transverse momentum of the produced 
gluon, k±<Q Sl one can neglect the logarithmic dependence on \x — y\ of the gluon distribution in Eq. (3) 
making analytical evaluation of the integral in Eq. (11) possible [17,2l]7 To calculate the integrals in Eq. (11) 
in this logarithmic approximation we need to estimate the following integral 

Hk, 2SV)= I d 2 z e-^ j^-^ e"* ^^^^ ^^ . (12) 

J \z-x\ 2 

After redefinition of variables z — > z + x the integral of Eq. (12) could be rewritten as 

/•OO 7 />27T 

I(k,x;y) = ie-^-5^) 2 ^V fe / — e -i* 2 ^/ d(f)e -ikzco S 4>-iz\x-y\Q^o S (<i>-0) (13) 



with <f> the angle between z_ and k and (3 the angle between z and x — y. The integration over <fi in Eq. (13) 
could be done (see 3.937.2 in [35]) yielding 



dz 



Hk,x;y) = 27r»e-**-*Ca-!0 Q ^ k / ™<rW I (zJ-X 2 ), (14) 

Jo z v 

where we have introduced a two dimensional complex vector 

A = k- t -Ql(x-y). (15) 

Differentiating Eq. (14) over k and integrating over z we finally obtain 

I(k,x;y) = -27T* A (i _ e -A 2 /Q 2 ) e -**-§ (*-v) 2 Ql. ( i 6) 

Inserting Eq. (9) into Eq. (11) and employing Eq. (16) yields 
d^diff a s C F 
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cPkdy 
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— / ^— ^fe-x £>~^h'y 
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where A* is the complex conjugate of A. 

Eq. (17) together with Eq. (15) provides us with the expression for the diffractive gluon production cross 
section for quarkonium-nucleus collisions explicitly in terms of gluon's transverse momentum k and the 
transverse separation of the onium state x — y. After integration over angles between k and x — y Eq. (17) 
could be easily rewritten in terms of the invariant mass of the produced particles which is approximately 
equal to M\ w fc 2 /x where y = In 1/x. 

To the degree that a quarkonium state can serve as a model of a baryon projectile such as proton, Eqs. 
(11) and (17) could describe the diffractive gluon production cross section in pA collisions. However, a much 
better model of pA collisions would involve a realistic proton consisting of three valence quarks and we are 
going to address this problem in Sect. IIB. 



B. Proton— Nucleus Collisions 

Diffractive gluon production in proton-nucleus collisions is different from our previous discussion of onium- 
nucleus scattering only by the fact that the proton consists of three valence quarks instead of quark-antiquark 
pair. Again we assume that the three quarks in the proton are in a color singlet state by themselves, i.e., 
there is no soft non-perturbative gluons in the proton's wave function to modify this picture. 

Similarly to the quarkonium case diffractive gluon production cross section could be expressed as a con- 
volution of the proton's wave function and the diffractive cross section of a qqq state on the target nucleus: 
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(18) 



with x 13 ,x 23 the distances between valence quarks in the proton and Qi,a2 light cone fractions of the 
proton's momentum carried by two of the quarks. 

The diagrams describing the diffractive gluon production in pA are shown in Fig. 3. They are similar to 
the onium-nucleus scattering diagrams of Fig. 1 with three quarks instead of qq pair in the projectile. The 
diagrams in Fig. 3 demonstrate that in order to calculate diffractive gluon production cross section we need 
to know the diffractive "propagators" of the qqq and qqqG states through the nucleus. 
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FIG. 3. Diffractive gluon production in proton-nucleus scattering (pA). 

First we are going to calculate the "propagator" of the qqq state through the nucleus, which is needed for 
the diagram in Fig. 3B. Similar to what was done in Sect. II A we need to consider scattering of three quarks 
on a single nucleon. The contributing diagrams are shown in Fig. 4. There are two types of diagrams in 
proton-nucleon interaction: there are diagrams where both gluons connect to the same quark line, as shown 
in Fig. 4A and there are diagrams where gluons connect to different quark lines, as shown in Fig. 4B. The 
color singlet configuration of three quarks is achieved by antisymmetrization over their color indices with 
the help of e a/ g 7 . 
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FIG. 4. Diagrams contributing to the interaction of the proton with a nucleon in the target nucleus. 



The diagram in Fig. 4A easily yields 



A 



-m^- 



(19) 



The diagram in Fig. 4B is a little harder to calculate. The color factor there is given by 



— e a/ 3>y {T a )f]p, (T a ) 7 y . 



Using the Fierz identity 
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we readily evaluate the diagram in Fig. 4B to be 
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(20) 



(21) 



To obtain the proton's propagator out of Eqs. (19) and (21) we have to sum over gluons' connections to all 
valence quarks in the proton and multiply the result by 2\fW~- ~W pa. The final expression for the proton's 
propagator is 



N e 



exp 



4 (27V C ) 



1 Nc 

\2 Z_^ '— « ~~ —3> 



2/^2 



»j= i;»^i 



which for the case of three colors becomes 



exp (~ [fe - x 2 ) 2 Q 2 s + {x 2 - x 3 ) 2 Q 2 s + (x, - x 3 ) 2 Q 2 s ] J 



(22) 



(23) 



Calculation of the qqqG "propagator" in the nucleus is more complicated because, unlike the above cases, 
the scattering of a single nucleon is not diagonal in the color space. Thus instead of exponentiating a scalar 
we will have to exponentiate a matrix in the color space. (For similar calculations of different quantities see 
[31].) To determine the color matrix we have to again consider the scattering of the qqqG system on a single 
nucleon. First let us note that the qqqG system can arise due to gluon emissions off each of the valence 
quarks, giving rise to a different color structure in each case. The possibilities are outlined in Fig. 5. Quarks 
there are labeled 1, 2, 3 carrying color indices a, (3, 7 correspondingly. 
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FIG. 5. Three possible ways of generating a qqqG fluctuation in a proton. 

The three states depicted in Fig. 5 correspond to three color matrices, which we will label v\, V2 and 113 
such that 



v\ = —j=e a >t3t(T a )aa>, 
V6 



v 2 = —7= e Q/ 3'7 (T a )pi3>, 



v 3 = -j=e a ^' (T a ) 7 y. 



(24) 



The factor of \j \f§ is arises in Eq. (24) due to (e Q/ 3 7 ) 2 = N c l — 6 and is included to factor out the color 
structure of the proton's wave function. 



As we will see below interactions with nucleons only transform each of the states in Eq. (24) into linear 
combinations of the others. It is therefore convenient to make an orthonormal basis out of these matrices. 
We note that since due to SU(N C ) group properties 



vi + v 2 + v 3 = 



(25) 



the matrices of Eq. (24) are not linearly independent, 
going to choose the following orthonormal basis 



1 



U-2 = V 2 + ~ V3, 



(In fact they represent the roots of SU(N C ).) We are 



"3 



V3 



v s , 



(26) 



where we have explicitly inserted N c = 3 and Cf = 4/3. Since in Fig. 5 we explicitly consider a baryon 
(proton) consisting of three valence quarks our results from here to the end of the section will be derived for 
N c = 3. To generalize this discussion to an arbitrary N c one would have to introduce iV c different matrices 
v±, . . . , i>at c and work in the N c — 1 dimensional space of the orthogonal basis u 2l ■ ■ ■ , un c - Doing that is 
beyond the scope of this paper. 
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FIG. 6. Interaction of the qqqG state with a nucleon. 

Let us now consider interaction of the qqqG state with a single nucleon in the nucleus. The corresponding 
diagrams are depicted in Fig. 6. Each of the diagrams is a matrix in the u 2 , u% linear space. To calculate the 
diagrams of Fig. 6 we have to first calculate the action of each of them on v 2 and ^3 and then use Eq. (26) 
to rewrite their contributions in the (u 2 , U3) basis. The result of a rather tedious calculation yields 
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(27b) 
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(27c) 



".-J*ta-*>($| Vf) (27e> 

M = ii> fe - a ,(_^-Vf) 

% = ^-^("j/ 9 2/ 9 ) (27h) 

where in Eq. (27a) we summed over gluon connections to all three valence quarks. Summing up the contri- 
butions of Eq. (27) and multiplying them by 2\/R 2 — b 2 pa yields 



-M(z) = 2\/R 2 -b 2 pa(M A + ... + M H ) (28) 

with the matrix M(z) given by 

(\Cz + T2 (C2 + Cl) + M (X23 + X13) - g X12 ^ (-C2 + Cl + X23 - Xis) \ 

Uj ^ ^3(-C2 + Cl+X23-Xl3) | C3 + i(C2 + Cl)-^(X23 + Xl3) + 1X12^1 [ ' 

where we have defined 

Q^\(z-x t fQl (30) 

and 

Xij = \{^-^?Ql (31) 

The "propagator" of the qqqG state through the nucleus is given by 

cxp[-M(z)], (32) 

which is also a matrix. 

Now we are in a position to write down a compact expression for the diffractive gluon production cross 
section in the proton- nucleus collisions. The initial proton state in the amplitude and in the complex 
conjugate amplitude could be in either one of the three states, v\, V2 and V3, which are given by 

vf = (-1 , ~), v T 2 = (1 , --L), vl = (0 , A) (33) 

in the (112, W3) basis. If the proton in the amplitude is in the state Vi and in the complex conjugate amplitude 
it is in the state Vj we can define a diffractive overlap functions 

11- = v f f e - M i£i) _ e -(2/9)(xi2+Xl3+X23) N ) ( e -M(z 2 ) _ e -(2/9)( Xl2 +Xl3+X23) N ) v .^ (34) 

where, as in Sect. IIA , z x and z_ 2 are the gluon's transverse coordinates in the amplitude and in the complex 
conjugate amplitude. The second exponent in Eq. (34) comes from Fig. 3B with the propagator of the proton 

10 



in the nucleus given by Eq. (23), which is a unit matrix in the (w2>W3) color space. The diffractive gluon 
production cross section in pA can be written in terms of IIjj's as 

j~qqqA 3 3 

da diff _ a s f 2 2 , 2 „ v^v^ Zi-SSi Z-2-Z.j 



*b* zi ** y: e ^-7]i • ^-?v n - ^ 



cPkdy (2^ tt „ 

Eqs. (29), (33), (34), (35), together with Eq. (18) provide us the answer for diffractive gluon production 
cross section in pA collisions. The matrix expression for IL,j in Eq. (34) could be rewritten in the usual 
"scalar" form. After some lengthy algebra one arrives at 



4 
n33 =3 
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^ UlJ V3 x/£(U 2 ; 
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V3 v U2; ^ ; ^ y^y ■ 2 

_ _L e -f (X12+X23+X13) + _L e -i [Mn^HM^)] fi( A j (3?) 

v3 V3 / 

where we have introduced 

D(z) = 4 [M 12 (z)] 2 + [M n (z) - M 22 (z)] 2 , (38) 



R( Z ) , cosh « + ^nfe)-M 22 fe) smh VM , (39 ) 

2 /D(i) 2 v ; 

L(*) = cosh ^M - M "^ 2 ^ sinh JM, (40) 

2 /D(iJ 2 v ; 

and Mn(z), M 22 (z) and Afi 2 (z) are the components of (symmetric) matrix M(z) given in Eq. (29). All 
other diagonal (Uu) and off-diagonal (Ily, i ^ j) propagators could be obtained from Eqs. (36) and (37) 
correspondingly by appropriate permutations of the indices of £'s & n d x' s - 

We have calculated above the diffractive gluon production cross section for quarkonium-nucleus and proton- 
nucleus collisions in the quasi-classical approximation (Eqs. (11) and (35)). The calculation for the case of 
pA collisions is much more sophisticated than the calculation for onium-nucleus scattering. This is due to the 
fact that the object in question is the diffractive cross section and therefore interactions with spectator quarks 
have to be included. When one calculates the total inclusive gluon production cross section the interactions 
of the target nucleus with the spectator quarks (or gluons) in the projectile disappear due to real-virtual 
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cancellations [17-19,31]. All one has to do there to obtain the answer is to calculate the interaction of a 
single projectile quark with the nucleus and then convolute it with the quark distribution function in the 
proton [17-19, 31]^. This is not the case for diffractive scattering: since the final state here is not "inclusive", 
i.e., we imposed a rapidity gap condition on the final state, only virtual exchanges contribute to interaction 
with the target, as seen in Figs. 1 and 3. Therefore there is no real-virtual cancellation for interactions with 
spectator quarks (quarks off which the gluon was not emitted neither in the amplitude nor in the complex 
conjugate amplitude). Interactions with all partons in the proton's or onium's wave function have to be 
included. That way, since our results of Eqs. (11) and (35) were derived for the models of hadrons consisting 
of valence quarks one should use them to describe the experimental data with caution. Realistic hadron's 
wave function at high energy may include non-perturbative fluctuations producing the so-called "intrinsic" 
sea quarks advocated by Brodsky et al in [36] and quantum fluctuations due to perturbativc QCD evolution 
producing ("extrinsic") sea quarks and gluons. Inclusion of either one of the two effects would tremendously 
complicate the diffractive gluon production cross section calculations presented above. Inclusion of the 
effects of sea quarks and gluons in diffractive production cross section is an important question which has 
to be addresses in a separate study. At the moment we may argue that perturbatively generated sea quarks 
and gluons would bring in higher powers of the strong coupling constant a s and could be neglected in the 
quasi-classical approach employed here. 

III. DEEP INELASTIC SCATTERING 

We can generalize the results of Sect. II A to the case of diffractive gluon production in the deep inelastic 
scattering (DIS) . The diagrams important for the small- a; gluon production are the same diagrams as shown 
in Fig. 1 only with the wave function of a virtual photon splitting into qq pair instead of onium wave function. 
The expression for the diffractive gluon production cross section in DIS is 

j -y'A^qqGA dn q ' iA 

^I_._L/ #ri ^« fc =,^ fc) , ,41, 

where da^fff/cPk dy is given by Eq. (11), r = x — y is the transverse separation of the quark- antiquark pair 
and z is the virtual photon's light cone momentum fraction carried by the quark. The square of the virtual 
photon's wave function $ 7 ^ qq (r, z) is a well known function and could be found for instance in [30]. 

An interesting question to address is the calculation of the inclusive gluon production cross section in DIS. 
This corresponds to the gluon production cross section without any restrictions on the final state of the 
target proton or nucleus. The problem is somewhat different from the gluon production in pA as discussed 
in [17,18,31]. In the model of the proton considered in [17,18,31] the interacting valence quark in the proton 
could have an arbitrary color in the initial state due to the soft intrinsic partons which would randomize the 
colors of valence quarks. Thus an independent summation was performed over the colors of this interacting 
quark. This could not be done in DIS: due to perturbative nature of the qq state at the leading order in 
a a there is no quantum fluctuations in the qq wave function randomizing the colors of the quark and the 
antiquark. Instead, in quasi-classical approximation employed here the quark and antiquark are in the color 
singlet state and an independent summation over the colors of one of them is not possible. Therefore we 
have to consider the diagrams with the gluon emitted off both the quark and the antiquark. 

There are two types of diagrams contributing to inclusive gluon production cross section in DIS. The 
diagrams may have the gluon emitted off the quark (antiquark) both in the amplitude and in the complex 
conjugate amplitude (symmetric diagrams) or they may have the gluon emitted off the quark (antiquark) 
in the amplitude and off the antiquark (quark) in the complex conjugate amplitude (asymmetric diagrams). 
The gluon production contribution in the diagrams of the first kind is exactly the same as in pA [17] (for a 



^This statement depends on the model of the proton: if we assume as above that the three valence quarks are 
initially in the color singlet state then the produced gluon could be emitted off different quarks in the amplitude and 
in the complex conjugate amplitude (see Sect. III). Nevertheless interactions with the quarks off which the gluon 
was not emitted cancel. 
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more detailed review see [21,31]). The contribution of the second (asymmetric) type of diagrams differs from 
the first in: (i) the interference terms, where the gluon is emitted before the collision in the amplitude and 
after the collision in the complex conjugate amplitude and vice versa; (ii) the late emission time terms, where 
the gluon is emitted after the interaction both in the amplitude and in the complex conjugate amplitude. 
To understand (i) let us discuss a specific case of a gluon emitted off the quark line before the collision in 
the amplitude and absorbed by the antiquark line after the collision in the complex conjugate amplitude. 
There the interactions with the target nucleons in which at least one of the exchanged Coulomb gluon lines 
connects to the quark line cancel by real-virtual cancellation with the diagrams where that gluon line is on 
the other side of the cut. Only the interactions with the gluon and the antiquark lines survive and could 
be easily resummed along the lines outlined above and in [17]. The late emission time diagrams (ii) in the 
case of asymmetric gluon emission are different from their symmetric counterparts because the real-virtual 
cancellation of interactions here does not happen. The reason is the color factors which are now different for 
the "real" and "virtual" diagrams. Similar to [17] one has to consider the last nucleon interacting with the qq 
pair before the gluon emission. With the help of Fierz identity the interaction diagrams for this last nucleon 
could be resummed and prove to be diagonal in the color space, which simplifies their exponentiation. At 
the end the one gluon inclusive production cross section for DIS could be written as 



, -y*A->qqGA 

Ua incl 

dPkdy 



with 



Hh-M> 



(XsCf 1 



(27T) 



1 

2^2 



d rdzQ 1 



da qqA 



d 2 bd 2 z 1 d 2 z 2 e-^^ 1 -z 2 ) 



d 2 kdy 



U*i 



(42) 



-(z 1 -x) 2 Q 2 /4 



-(2 2 -x) 2 Q^/4 + e -(2 1 - 22 ) 2 Q=/4 N j + 



ii -y 



ki - y\ 2 \z 2 - y\ 



- 2 ~4d 



-( 2l -y) 2 Q 2 /4 _ e -(£ 2 -y) 2 Q 2 /4 



+ e -(£ 1 "£ 2 ) 2 Q 2 /4A _ Jl 

J \z, 



z-y-x z 2 -y 



\z.i - x\ 2 \z 2 



( e -(2L-y) 2 Q 2 Ji _ e -(z 1 -y) 2 Q 2 s /4 _ e -(z 2 -s.) 2 Q 2 j4 + e -(^^ 2 ) 2 Q 2 s/A 



s.i-y 



\z x -y\ 2 \z 2 -x\ 2 



e -(x-y) 2 Q 2 j4 _ e -( 2l -x) 2 Q 2 /4 _ e -(z 2 ~y) 2 Q 2 /4 + £ - (^- z 2 ) 2 Q 2 / 4 



(43) 



As one can easily check the cross section of Eq. (43) goes to zero in the limit of zero dipolc size x — > y. 

Eqs. (42) and (43) provide us with the one gluon inclusive production cross section in DIS calculated 
in the quasi-classical approximation and could be used to describe the pion or minijet production data at 
HERA. 
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